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Abstract 

We introduce a theoretical approach based on scattering theory and total 
energy methods that treats transport non-linearities, conformational changes 
and charging effects in molecular wires in a unified way. We apply this ap- 
proach to molecular wires consisting of chain molecules with different elec- 
tronic and structural properties bonded to metal contacts. We show that 
£h ' non-linear transport in all of these systems can be understood in terms of a 

single physical mechanism and predict that negative differential resistance at 
high bias should be a generic property of such molecular wires. 
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I. INTRODUCTION 

A molecular wire in its simplest definition is a single molecule that forms an electrically 
conducting bridge between a pair of metallic contacts. During the last few years molecular 
wires have been realized in the laboratory and the electric current that flows through an 
individual molecule in response to a potential bias applied between the contacts has been 
measuredaQ. The current-voltage characteristics of molecular wires are sensitive to the 
electronic structure of the molecule and contacts, which depends on the spatial arrangement 
(conformation) of the atoms that make up the wire and the distribution of electronic charge 
along the wire: Molecular wires have been observed to switch between conducting and 
non-conducting states due to a voltage-induced redox reaction, i.e., a change of the charge 
distribution in the wireB, and also in response to bias-driven conformational changesa. These 
experiments indicate that transport non-linearities, conformational changes and charging 
in molecular wires are mutually interdependent phenomena, each of which can strongly 
influence the others, and that they should therefore be treated theoretically in a unified 
way. Theoretical work to date has studied some of these phenomena independentlycll3ll2l 
but there has been little work directed at bringing all three of them together within a single 
theoretical framework. 

In this article we introduce a tractable theoretical approach that unifies the three phe- 
nomena and present the results of self-consistent calculations for the simplest molecular 
wires, namely, small chain molecules bonded to metal electrodes. At low bias these molecu- 
lar wires are conductors or insulators depending on which atoms make up the chain molecule 
and electrodes and on the number of atoms in the chain. Remarkably, we find that the non- 
linear transport characteristics of these apparently dissimilar systems are strikingly similar 
and can be understood in terms of a single physical picture: We find that in the non-linear 
transport regime the molecular wire acquires excess electronic charge. This charge resides 
in a molecular orbital which, being partly occupied, lies at the Fermi level and therefore 
controls electron transport through the molecule. (If there is no molecular orbital near the 
Fermi level at zero bias, we find that the conformation of the molecule changes at finite 
bias in such a way that a suitable orbital appears there). As the bias increases, the partly 
occupied orbital acquires more charge and the electrostatic contribution to its energy rises 
together with the electro-chemical potential of the source electrode. When this orbital rises 
well above the electro-chemical potential of the drain electrode the electron flux that it 
transmits from drain to source approaches zero and ceases to decline with increasing bias 
while the flux that it transmits from source to drain cannot increase further significantly. 
Thus the net current through the molecular wire saturates. At the same time a charge 
density wave grows along the wire and atomic displacements occur, lengthening and thus 
weakening certain molecular bonds. We find the bond weakening to be sufficient that after 
saturating the current carried by the wire decreases as the bias increases further. We predict 
this novel negative differential resistance phenomenon to be a generic property of molecular 
wires that involve chain molecules. 

Sec. II presents the theoretical formulation of the model that includes charge transfer 
and conformational change in the non-linear transport regime. The transport characteristics 
of two simple linear chain-like molecules (homogeneous and inhomogeneous) are presented 
in Sec. Ill based on calculations using the model of Sec. II. Finally, Sec. IV summarizes the 



results of this work. 



II. THEORY 

Recent theories of electronic transport in molecular wires have been based on semi- 
empirical tight-binding models (such as the extended Huckel model) and on ab-initio com- 
putations that are founded on density functional theory. For standard values of the tight 
binding model parameters adopted from the quantum chemistry literature, the tight bind- 
ing models have yielded results that are similar to (and often in semi- quantitative agree- 
ment with) those obtained from the ab-initio computations. For example, tight bind- 
ing calculations^-!] have yielded current-voltage characteristics for the benzene-dithiolate 
molecule bonded between gold nano-contacts that are similar to those that have been found 
in ab-initio calculations! These studiesQ'tll have made it evident that the differences be- 
tween the results of the tight-binding and ab-initio calculations for this system are less 
significant than the uncertainties that are due to the present limited understanding of the 
conformational aspects of bonding between the molecule and metal contacts. Tight binding 
calculations of electronic transport through molecular wires consisting of chains of carbon 
atoms connecting metal contacts have also yielded resultscJ that that agree remarkably well 
with those found in ab-initio calculations for the same systems^. The ab-initio calculations 
have the advantage of not relying on the use of empirical parameters. However they are 
much more computationally intensive than tight-binding calculations. Therefore they are 
limited in their application to smaller molecular systems and rely on the use of the jellium 
approximation to treat the metallic contacts. For the same reason it is difficult to study 
conformational changes using the ab-initio techniques and the molecules were assumed to be 
rigid in the above ab-initio studies.Q'tj Since the ab-initio calculations are based on density 
functional theory which is known to be rigorously valid only for ground state properties, 
their validity for transport calculations far from equilibrium is uncertain a priori. Further- 
more even in calculations of ground state properties they rely in practice on the use of the 
local density approximation. However, it is reasonable to expect them to be more reliable 
than calculations based on simpler mean field theories such as the Hartree-Fock approxima- 
tion. While the Keldysh technique can in principle provide a more rigorous approach to the 
problem of transport far from equilibrium, its practical application to real molecular wires 
still presents major challenges at this time. 

A shortcoming of the tight-binding approach to molecular wire transport has been that it 
has not provided a way to calculate the conformational changes that occur in a molecular wire 
when an electric current is driven through it, or to treat the influence of these conformational 
changes and of charging of the molecule on the current in a consistent way, within the same 
theoretical framework. The formalism that we present below accomplishes this, however 
it should be emphasized that it, like the ab-initio calculations mentioned above, is based 
on density functional theory and should therefore also be regarded as a practical mean 
field approach to the difficult many-body problem of transport in molecular wires far from 
equilibrium, and not as a rigorous treatment of this problem. 

In our tight-binding formalism we represent the electronic energy S of the molecular wire 
as a functional £[{iVj}] of the of the numbers iVj of electrons occupying atomic sites i at 



positions R^. 

8[{Ni}] = £[{iVj] -eJ2N l V ext (R t ) + 

2 4vreo^|^-^.| [1 * tJ 

This is a tight-binding analog of the Kohn-Sham!^! density functional: /C[{iVj}] is the tight- 
binding energy of a non-interacting electron system with the same electron distribution {Ni} 
as the interacting system and V ex t is the corresponding external potential. The next term 
is the electrostatic interaction energy between electrons on different sites. T represents the 
exchange- correlation energy and on-site interactions. In the spirit of the Kohn-Sham local 
density approximationEJ we approximate T in terms of a local function fi(Ni) that depends 
on the atomic species at site i as jFfjiVj}] = J2ifi(Ni). Varying £[{iVj}] with respect to 
{N} with the conformation {Ri} held fixed we obtain the tight binding analog 

Y, *r< + (<*» + ^im = e x ^l (2) 

jn 

of the Kohn-Sham equations for the effective one-electron eigenstates \1/ A = ifj^ m \i,m) and 
energy eigenvalues E x . t™* 1 is the hopping matrix element between orbitals m and n on 
atoms i and j, e im is the bare orbital energy, 

V ext includes potentials due to the bias applied to the contacts and image charges. 
5fi(Ni)/SNi is the variation of the electronic site energy with the electron occupation iVj 
of the site. In the numerical work reported below we approximate it by Sfi(Ni)/8Ni = 
(Ni — N^)Ui where N® is the electron occupation of the neutral atom and Ui is the on site 
interaction energy whose value we take from the semi-empirical chemistry literature where it 
is estimated from experimenttH. It is possible to also include corrections to this term that are 
of higher order in (iVj — N®) since the relevant empirical coefficients are available;E-3 we did 
not do so in our numerical work because they are small for the species that we considered. 
For a given molecular conformation {Ri} and applied bias we solve Eq.(^) self- 
consistently with Ni in Eq.(|3]) obtained by integrating J^ n \ipf n \ 2 over the occupied eigen- 
states of Eq.(0). The electric current carried by the wire is then given by Landauer-Biittiker 
theory: 

I(V) = j J dE T(E, V) [F(E, ti a ) - F(E, /i d )] (4) 

Eq.(f|) relates the current I at finite bias voltage V to the transmission probability T(E, V) 
for an electron to scatter from the source metallic contact through the molecule and into the 
drain at an energy E. F(E,n) is the Fermi function, and // Sj d = E F ± eV/2 are the source 
and drain electro-chemical potentials with Ep the common Fermi energy of the contacts. 



We obtain T(E, V) from the eigensolutions \1/ A of Eq.@ and the velocities of electron Bloch 
states in the contacts in the standard wayM. 

In the Landauer picture of transport that gives rise to Eq.(|), the electrons populate 
scattering states that flow from deep in the contacts towards the wire according to the 
Fermi functions F(E,fi a ) of the contacts. These states are partly transmitted through 
the wire giving rise to the current I and partly reflected. We adopt the same physical 
picture in calculating the actual conformation that the molecular wire takes in the presence 
of the applied bias and transmitted current: We assume that the asymptotic regions of 
the contacts constitute ideal leads where the electrons may be treated as non-interacting. 
We also assume that the filling of electron states that are incident on the wire from each 
contact is described completely by the Fermi function of that contact, i.e., it depends only 
on the applied bias and the temperature and not on the conformation of the wire. [In the 
present work the single particle states of Landauer theory are replaced by the Kohn-Sham- 
like effective single-particle states that are the self-consistent solutions of Eq. (0).] Then, if 
electron scattering at the wire is elastic, the conformation that the wire takes is that which 
minimizes the total energy of the system, the filling of the asymptotic incident states with 
electrons being held fixed. The total energy E to t{{Ri}) for any trial conformation {Ri} is 
found from the self-consistent solution of Eq. (g). Its electronic part E e is equal to £[{iVj}] 
for the conformation {Ri} (and the fixed incident electron population) and is given by 



E e = J E[D S (E)F(E, fi a ) + D d (E)F(E, fi d )]dE + 



e 2 v^ N,N. 
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where Ni is the self-consistent value of Ni that corresponds to the solution of (0) and D S}d are 
density of states factors for the asymptotic incident states in the contacts. Equation (|5J) is the 
non-equilibrium tight-binding analog of Equation (2.10) of Kohn and ShamEJ: The integrand 
in Equation (Bl) corresponds to the sum of the energy eigenvalues of the occupied Kohn-Sham 
effective single-particle states. These single particle states correspond to electrons that arrive 
at the molecular wire from the source and drain leads and are populated according to the 
Fermi distributions of the source and drain leads respectively. The first summation corrects 
for the double counting of the Coulomb energy due to interactions between electrons that 
is implicit in the integral that precedes it. The second summation corrects the contribution 
of the on-site energy to the total electronic energy. 

E e includes implicitly the interaction energy between the nuclei and the electrons of the 
wire but omits the electrostatic energy of the nuclei of charge eZi interacting with each other 
and with the potential V ext due to the applied bias and image charges in the contacts. Thus 
to find the actual conformation of the wire we minimize E tot ({Ri}) = E e ({Ri}) + E n ({Ri}) 
w.r.t. {Ri} where 

Notice that the energy eigenvalue integral in Eq. (||) is independent of the conformation 
of the wire. This independence is a consequence of working with an open system with a 



continuous spectrum and of our adoption of the Landauer picture of the contacts discussed 
above. 



III. PREDICTIONS FOR SOME SIMPLE MOLECULAR WIRE SYSTEMS 

We have applied this methodology to two classes of simple molecular wires consisting 
of chain molecules bonded to metal contacts (see Fig. [I]). The contacts are modeled as 
single-channel ideal leads with Fermi Energy Ep = —10 eV and the energy band ranging 
from —17 eV to —3 eV. They are assumed to form planar equipotential surfaces. For one 
class (inhomogeneous wires) the molecules are representative of finite polyacetylene chains; 
we use the model parameters of Su et alx^. Each atom of the backbone of the chain has a 
single orbital with bare energy ec = —11.4 eV, nearest-neighbor hopping t° c = —2.5 eV for 
the undimerized case and charging energy Uc = —11.5 eV. Each end of the chain is each 
bonded at a perpendicular distance of 2A over a single metal atom of a gold electrode 
which has a bare orbital energy 6m = —10 eV, a charging energy Um = 8.5 eV, and 
is coupled to the end of the chain with a hopping energy t^ M = —2 eV. For the other 
class (homogeneous wires) which is representative of metal atomic chains bonded to metal 
electrodes, the unrelaxed geometry is a linear chain with a 2.5A spacing. The same bare 
orbital and charging energies are used as for the above gold metal electrode atom, with a 
hopping energy t° M = —3 eV. As the atomic positions relax in the course of minimizing the 
energy, the hopping parameters t change. Following Su et al.tB we assume t — t° + a5X 
where 5X is the change in bond length. For C-C, M-C and M-M bonds, ac,c — 4.125, 
&c,m = 2.5 and olm,m = 3.0 eV/A. We calculate the current- volt age characteristics of 6 and 
7 atom inhomogeneous and homogeneous chains, allowing the atoms to relax both in the 
presence and absence of a bias voltage. 

In Fig. |2|, the calculated site charges — e(iVj — N®) are shown for bias voltages of 0, 1 and 
2V. For the homogeneous chains (Fig. |2|a,b) each atom is neutral at 0V. At higher bias, a 
charge density wave (CDW) forms in both the 6 and 7 atom chains. The CDW arises from 
a superposition of standing electron wavesc2lfill in the wire. The relaxed atomic positions 
are the lateral positions of the symbols in each graph. The vertical lines show the relaxed 
positions at 0V; the unrelaxed positions match the ticks on the z axis. At 0V the 6 atom 
chain dimerizes while the 7 atom chain has a frustrated dimerization pattern. At higher V, 
the weakening of some bonds amplifies the standing waves and CDWs. The inhomogeneous 
chains (Fig. 0c, d) show similar features although the CDWs are present even at 0V and 
there is a net transfer of charge from the metal contacts to these wires. The atoms at the 
ends of the chains in Fig. 0c, d are metal electrode atoms. These metal atoms are positively 
charged as charge is transferred from them to the chain. For the 6 atom chain, this charge 
transfer at the ends of the chain frustrates the CDW, and a charged soliton forms. At V 
the soliton is at the center of the chain. It moves away from the more negatively charged 
end of the chain at larger bias. The soliton is associated with a stretched bond. For both 
inhomogeneous chains, there is considerable weakening of the bond between the chain and 
drain at higher bias. This will be seen to have an important effect on the current flowing 
through the wire. 

These charging and conformational changes are intimately linked to the transport of elec- 



trons through the wires. The voltage-dependent transmission probabilities for 0, 1 and 2V 
are shown in Fig. |3|. Consider first the case of zero bias and odd atom chains (Fig. |^b,d). In 
isolation these chains have an odd number of valence electrons (in this case 7) which occupy 
molecular orbitals, the highest occupied molecular orbital (HOMO) being half-filled. In an 
open system if this orbital were well below Ep it would have to be completely filled and if it 
were well above the Ep it would be completely empty. Both situations lead to unfavorable 
charging of the wire. Thus the HOMO must be at Ep and transmission resonances due to 
the HOMO can be seen at Ep = —10 eV in the two graphs. (The HOMO of the inhomoge- 
neous wire is more than half full due to charge transfer from the contacts. Thus the center 
of the 0V HOMO resonance in Fig. ^d is slightly below Ep). The 6 atom inhomogeneous 
chain also has a resonance near Ep at zero bias (Fig. ||c). This can be explained similarly 
because there is charge transfer which must go to a partly filled state near Ep. This partly 
filled state leads to the charged soliton in Fig. 0c. However, the 6 atom homogeneous chain 
has no resonance at Ep at 0V (Fig. |3]a) since there is no unpaired electron and no charge 
transfer at zero bias. Indeed it is energetically favorable for this chain to dimerize at 0V 
pushing its HOMO lower in energy and widening the gap around Ep. This explains why 
there are resonances at Ep at zero bias for some of the systems but not for others. Charge 
transfer also explains why there are resonances at Ep for all of the chains at non-zero bias: 
All of the wires become partly charged at non-zero bias. Thus all of them must have a partly 
filled level at non-zero bias and this level must be near Ep. For the 6 atom homogeneous 
chain, a significant rearrangement of its atomic positions is required to generate this state 
near Ep. (The appearance of fewer resonances than atoms in the chain in some plots is due 
to degeneracy of some of the transmitting states). 

The current vs. voltage characteristics calculated from these transmission probabilities 
are shown in Fig. |j. Their physical meaning is as follows: At moderate and high bias the 
resonance near Ep carries the net current through the wire. As the bias increases the wire 
becomes charged and the Coulomb potential due to that charge raises the energy of the 
resonant state until the bottom edge of the resonance is above the electro-chemical potential 
of the drain. Before this occurs the resonance transmits electrons coming from both the 
source and drain. As the bias increases there are fewer electrons from the drain and more 
from the source going through the resonance so the current is growing. Eventually there are 
no electrons from the drain going through the resonance and the current saturates. It then 
declines and negative differential resistance sets in because the stretching and weakening of 
bonds inhibits electron transmission through the wire. 

IV. SUMMARY 

In summary, we have presented theoretical work that for the first time unifies non- 
linear transport, charging and conformational change in molecular wires. We have shown 
that the effects of charging and conformational change can be understood together in terms 
of a remarkably robust physical mechanism of non-linear transport and predict negative 
differential resistance at high bias. One class of molecular wires that we have modeled (the 
homogeneous wires) is representative of chains of metal atoms connecting metallic contacts 
such as have recently been realized in the laboratory for the case of gold.£3'E3 Experiments 



carried out to determine whether these systems exhibit negative differential resistance at 
high bias would be of interest. 

We thank R. Hill for rewarding discussions. This work was supported by NSERC and 
by the Canadian Institute for Advanced Research. 
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Emberly and Kirczenow Fig. 1 

FIG. 1. Schematic of atomic wires bonded to ideal metallic leads, (a) Atomic configuration of 
homogeneous wire, (b) Atomic configuration of inhomogeneous wire (note the inclusion of metal 
atoms on the ends to model charge transfer between the metallic contacts and the atomic wire). 
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Emberly and Kirczenow Fig. 2 

FIG. 2. Site charges —e{Ni{Ri) — Nf) in units of e and atomic positions (A) for different 
voltages. Circle = 0V, square = 1 V, triangle = 2V. (a), (b): 6 and 7 atom homogeneous wire. 
(c), (d): 6 and 7 atom inhomogeneous wire. 
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Emberly and Kirczenow Fig. 3 

FIG. 3. Voltage-dependent transmission probability at (solid line), 1 (dashed) and 2V 
(dash-dotted). (a),(b): 6 and 7 atom homogeneous wire, (c), (d): 6 and 7 atom inhomogeneous 
wire. 
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Emberly and Kirczenow Fig. 4 

FIG. 4. Calculated current vs. voltage for (a) homogeneous wires and (b) inhomogeneous wires 
at T = 100 K. 
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